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In this paper we study brane-world scenarios with a bulk scalar field, using a covariant formalism 
to obtain a 4D Einstein equation via projection onto the brane. We discuss, in detail, the effects of 
the bulk on the brane and how the scalar field contribute to the gravitational effects. We also discuss 
choice of conformal frame and show that the frame selected by the induced metric provides a natural 
choice. We demonstrate our formalism by applying it to cosmological scenarios of Randall-Sundrum 
and Hofava-Witten type models. Finally we consider the cosmology of models where the scalar 
field couples non-minimally to the matter on the brane. This gives rise to a novel scenario where 
, the universe expands from a finite scale factor with an initial period of accelerated expansion, thus 

• avoiding the singularity and fiatness problem of the standard big bang model. 
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■ I. INTRODUCTION 

' The standard Big-Bang model based on the FRW spacetime now has many observational successes, notably the 
1 observed expansion of distant galaxies, the synthesis of light elements during Big Bang Nucleosynthesis (BBN) and 
I ■ the isotropy of the Cosmic Microwave Background (CMB). These successes mainly probe what, from the point of 
view of this paper, we shall describe as the late universe, that is, after neutrino decoupling — which is just before 
BBN. Our knowledge of the universe before BBN is very limited and hence it is possible for substantial modifications 
to our understanding of the expansion of the universe during such an epoch. 

Motivated by attempts to solve the hierarchy problem of particle physics and to reconcile cosmology with string 
theory, much interest has been generated by the idea that this Big-Bang model is embedded inside some higher 
; dimensional spacetime; the 4D universe that we see being thought as a 3-brane or brane-world. In such models the 

■ matter fields of the standard model are localized by some mechanism to only have support on the 3-brane, with gravity 
, and other undetected particles, such as scalar fields, propagating in the directions perpendicular to the brane or 'bulk' 
' as they have become known. In the most of popular of these models, the Randall-Sundrum (RS) scenario 1^,^, the 
' brane is a hypersurface (or co-dimension 1 brane) in a 5D spacetime with Z2 symmetry along the extra dimension. 

' In fact, RS made two simple proposals (see also ref. ||^ for related earlier work): the first ||] (RSI) has two branes 
O I of opposite tension in an anti-de-Sitter (adS) background spacetime with Z2 symmetry. In such a model the global 
D spacetime structure is warped, that is, the spatial 3-metric contains an exponential or warp factor which is a function 
of the coordinate along the extra-dimension. They suggested that this could possibly solve the hierarchy problem of 
^ , particle physics by relating the Planck scale on our brane (the one with negative tension) to that on the other via 
an exponential function of the distance between them. Hence, a large hierarchy (^ 10^^) can be generated from a 
. , relatively small dimensionless number. 

• In the second such model Q (RS2), there is only a single positive tension brane: it being possible to consider this 
as a special case of RSI with physical interpretation of the branes reversed and the negative tension brane taken 
off to infinity. In this scenario, which has an effectively infinite extra-dimension cut-off by the adS horizon, they 
showed that there exists a normalizable gravitational zero mode coupled to the brane which implies that to some 
degree the gravitational effects on the brane are close to those predicted by Newtonian theory and Einstein's 4D 
General Relativity (GR), even though the spacetime has an extra dimension. More precisely they showed that the 
gravitational theory experienced by observers on the brane is a tensor theory, being mediated by a massless spin 2 
particle, with a tower of massive Kaluza-Klein type states due to the extra dimension. Although there is no natural 
solution to the hierarchy problem in this model, much interest has been generated by the possibility of creating such 
a model within a modern fundamental theory based on string or M-theory. 

Interestingly, a very similar proposal already existed in the literature: the Hofava and Witten (HW) theory 
In this model, which is derived from the strong coupling behaviour of E% x E% string theories in lOD, the Z2 symmetry 
is natural; it being a consequence of the orbifold structure of the related IID M-Theory. Orbifolds, manifolds which 
are smooth except for a finite number of points, occur naturally in realistic compactifications of string theories since 
compactification on smooth manifolds cannot generate certain aspects of the Standard Model, for example, the chiral 
properties of neutrinos. Furthermore, the fixed points of the orbifold give natural places for the branes to live. In the 
HW case the orbifold is /Z2 so there are two fixed points. This is similar to the RSI model except that vacuum 
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energy contributions are generated by the effects of a scalar field with bulk support. Simple cosmological solutions 
were studied, in refs. (8|J|], although very little was learned about the detailed dynamics of the expansion. 

The existence of the gravitational zero mode coupled to the brane suggests that the gravitational effects are close to 
those of 4D GR, but there are subtle differences, particularly on small and large scales. For example, using a linearized 
analysis, it was shown in ref. pO[ that a localized mass distribution on the positive tension brane of RSI will effect 
the gravity on the negative tension brane in such a way as to make it appear to be like 4D Brans-Dicke (BD) gravity. 
The fact that any brane-world model must give rise to a universe which expands in an essentially similar way to the 
standard FRW universe, at least from the epoch of BBN, implies that one can probe the gravitational dynamics of 
such model at large distances using cosmology. Probing the viability of brane-world models in which there is an extra 
scalar field with bulk support using cosmology is the main motivation of this paper, along with the possibility that 
such scenarios might alleviate the many naturalness problems, such as the cosmological constant problem, prevalent 
in 4D FRW models. 

A large number of papers have attempted to understand the cosmological expansion rate in RS type|^ brane-world 
models, this being the first thing that one must get right before making any more exotic predictions. Here, we give 
a brief, and no doubt incomplete, history of the main recent developments. One of the remarkable aspects of this 
area is that it has been shown that there exist spatially homogeneous and isotropic brane solutions whose expansion 
can be understood almost totally without knowledge of the detailed 5D solution, with only small corrections to the 
expansion rate due to the effects of the bulk. Initially, it was suggested that a simple model for a 3-brane in a 

flat 5D background would lead to a Hubble expansion rate H (x p, where p is the density of the matter on the brane, 
in complete conflict with the observed expansion rate {H^ oc p in a 4D FRW universe) and hence with BBN. However, 
this ignored the lesson of the RSI and RS2 models. In the simplest versions of these two models vacuum energy, or 
cosmological constant (CC), contributions were placed on the two branes along with a negative CC in the bulk to give 
asymptotic adS geometry. Decomposing the matter on the brane into a CC plus other matter and including a negative 
CC in the bulk rectified the problem assuming the standard RS relation between the CCs [p^pT[. This modification 
leads to the rath er perp lexing result that on the negative tension brane cx ~p. It was sugg ested in and 
more recently in poy2l[ that an averaging procedure a la Kaluza-Klein could rectify this situation if one introduced 
a bulk scalar field to stabilize the extra dimension via, for example, the Goldberger-Wise mechanism |2^. Another 
interestingconsequence of these models is a radiation type contribution to the expansion rate due to the effects of 
the bulk 

Cosmological solutions in brane-world models with a bulk scalar field (see refs. [p5|-[27|] for previous work on this 
subject) are necessarily more complicated due to the extra field. However, they are probably more interesting since it 
would seem unlikely that only gravity would propagate in the bulk. In this paper we attempt to make some progress 
in understanding how the inclusion of such a field will effect the cosmological evolution, and in particular whether the 
standard FRW expansion is possible at late times. In contrast to the case of an RS type brane-world model, we find 
that there is an obstacle to a closed form solution in these models since the evolution on the brane can be affected in 
an almost arbitrary way by the dynamics of the scalar field in the bulk. In particular, it is not possible to derive an 
effective equation for the evolution of the scalar field on the brane, in terms of quantities which are only defined on 
the brane. In order to make any statements about the late-time behaviour of the expansion, therefore we are forced 
to make the reasonable assumption that the value of the scalar field on the brane is either stabilized or slow varying. 
We will see that the expansion rate is related to value of the scalar field, and hence if it is varying, one would have 
a model with a variable effective gravitational constant, which is constrained by observations to be at least slowly 
varying after BBN. An interesting consequence of this approximation is that the cosmological expansion of stabilized 
or slowly varying dilatonic brane worlds is essentially that of an appropriately tuned RS type model. 

Our approach to this problem will be to use the Gauss-Codacci formalism which gives a coordinate independent 
description of the dynamics on the brane in terms of effective 4D Einstein equations [ pT| . In section || we extend the 



analysis of ref. 1 17| to the situation where there is a scalar field in the bulk. Note that these methods are very different 
to those used in most of the other literature [pT|- |T6| , p3|j2^ ] probing cosmology on RS type brane-worlds where specific 
form of the metric is used. We show how the results of these two different approaches can be reconciled easily and, 
as a consequence, gain an improved understanding of the effects of the bulk. As we are dealing with a theory with 
a scalar field, the possibility of choosing a different conformal frame arises, as in 4D scalar-tensor gravity theories. 



We discuss the possibility of conformal scalings of the metric in section [ID. We show, firstly, that our frame is 



the preferable one on the grounds that it is the only frame in which the energy-momentum contribution due to the 



'^We shall refer to a RS type model as one with just constant vacuum energy terms, as opposed to those which are generated 
by the existence of a scalar field. 
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scalar field does not contain second derivatives of the scalar field and, secondly, that the conformal scaling, although 
allowing the sign of the Hubble parameter to change, does not affect the relationship between the signs of the brane 
tension and the effective Newton constant. 

In section III we apply these methods, first, to the simplest RS type model obtaining the results of ref. jl^, before 
applying them to the more complicated models which include a stabilizing scalar field potential for the size of the 
extra dimension. We show that once the scalar field becomes stabilized the 4D Einstein equations become exactly 
those of the RS model and hence a negative tension brane leads to oc —p as in the case of constant potentials. We 
discuss in detail the conflict between the results obtained by this approach and those of other analyses, which have 
concluded that the standard FRW result should apply. This raises the important question of what to interpret as the 
effective 4D theory. 

We then study a simple HW type model with exponential potentials in the bulk and on the brane. Assuming that 
the matter on the brane is not coupled to the scalar field, and the scalar field on the brane is varying slowly, we 
show that the RS type cosmological evolution is still valid, albeit now with an effective gravitational constant which 
is varying with time. It is interesting to note that in such models energy conservation on the brane is modified by the 
variation of the scalar field giving rise to a novel scenario for the variation of the coupling constants. If one allows 
the matter on the brane to be coupled to the scalar field as well, then it is possible to achieve accelerated expansion 
on the brane leading to an interesting bulk driven inflation model even when the scalar field is stabilized. 

II. FORMALISM 
A. Geometrical description of the spacetime 

The basic method is to foliate the 5D space-time in the direction of the extra-dimension as done in ref. . This 
makes a 4-1-1 split of the Einstein equations in a similar way to the standard 3-1-1 split used in treatments of the 
Cauchy initial value problem in GR (see, for example, ref. j28j). This geometrical split allows us to understand the 
precise gravitational dynamics on the brane without any possible problems with unphysical gauge degrees of freedom. 

More precisely, we wish to study space-time^ including a smooth 4D time-like hypersurface (or co-dimension 1 
brane) which we can express, in terms of the coordinates x"" with (a = 0,1,2,3,4), by /(x") = 0, where / is a 
real-valued function. Hence, we can define a space-like normal vector daf and a unit vector Ua parallel to it by 

na^daf/idbfd'fy/^ (1) 

One would like to extend Ua to follow a geodesic congruence, that is, so that it satisfies the geodesic equation, 
n''Vf,n° = 0. Near the brane, this will define a smooth foliation, but it may not be possible to foliate the entire bulk 
in this way. Moreover, in cases where there is more than one brane, a geodesic congruence which is normal to one 
brane may not be normal to another. For example, a coordinate system which is Gaussian-normal with respect to one 
of the branes, that is, Ua = (0, 0, 0, 0, 1) on that brane, is not necessarily Gaussian-normal with respect to another |^0|. 
(See also Q for a discussion of the regions covered by Gaussian normal coordinates in Schwarzschild-adS space.) 

One can define the metric, or (first) fundamental tensor, of the foliation, hat, in terms of the normal vector to 
be its orthogonal complement, 

h\ = S\ - n'^nt, . (2) 

As well as being the higher-dimensional manifestation of the induced metric on the brane, this is a projection tensor 
which, along with the it orthogonal complement, allows one to project other tensors tangential on to the brane. As 
an example consider any vector, y , which can be written as 

- hW' + {mV')na . (3) 

Clearly, the first part is tangential to the brane and the other is perpendicular. We shall see in subsequent sections 
that this projection property is an important concept when to attempting to understand brane-world since it allows 
one to split the dynamics up into those tangential to the brane and those perpendicular. 



^Note that, throughout, we will use the same sign conventions as |2^,^: in particular, we use a ( — h -|- -|- -I-) signature for 
the 5D space-time metric. 
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The extrinsic curvature, or second fundamental tensor, is defined to be 



Kab = Va^b = ha^hb^VcTld , (4) 

which is symmetric by Frobenius' theorem and tangent to the hypersurface. It can also be written in terms of Lie 
derivatives of the metric with respect to the normal 

2i^ah = £n9ab = £ nhab , (5) 

and its Lie derivative is given by 

£nKab^KacK^'-Racbdn'n\ (6) 

where Rabcd is the Riemann curvature tensor of the background space-time. In a Gaussian-normal coordinate system, 
the extrinsic curvature is half of the derivative of the hypersurface metric, hab, with respect to the coordinate normal 
to the hypersurface. 

Within this setup, we can use the Gauss-Codacci formalism to compute effective 4D Einstein equations for the 
gravitational dynamics on the brane. We will denote all curvature tensors pertaining to the hypersurface metric, /lab, 
with a bar over the letter, whereas those pertaining to gab will have no bar. Similarly, V is the covariant derivative 
that preserves hab whereas V is the covariant derivative that preserves gab- The Gauss and Codacci equations are 
then 

Rabcd = ah^bh^ ch"^ dRjklm + Z/CaJc-ftTdJb , (7) 

ybK\ - VaA' = n^h\Rbc = n^h\Gbc ■ (8) 
Following the approach of ref. |1^], we use the decomposition of the 5D Riemann tensor 128] 



2/ \ 1 

Rabcd = -^\ga[cRd]b — 9b[cRd]a I — QR9a[c9d]b + Cabcd ■ 



(9) 



which we substitute into (^. Contracting appropriately gives the following expression for the 4D effective Einstein 
tensor, 



Gab = ^ !^Gcdh\h\ + (^Gcdn'n'' - ^G^ habj + KKab - Ka^K^c " ^ " K^'^Kcd) hab - Eab 



(10) 



where Eab is the electric part of the Weyl tensor with respect to Ua, given by 

Eab = CacbdU^n'' . (11) 

We should note that at this stage everything we have derived is a geometric identity and hence true on all of the 
4D hyperspaces in the foliation (assumed to be well-defined). In subsequent sections we will use junction conditions 
to compute Kab close to the brane and use this expression for the Einstein tensor to make a comparison with the 
predictions for cosmological expansion of the brane and those of 4D GR and other gravity theories 



B. The action and variation 



In order to impose some dynamics on the geometrical identities of the previous section we will derive equations of 
motion from a variational principle which includes the standard Einstein-Hilbert term albeit in 5D, a scalar field (j) 
which has support in the bulk as well as on the brane, and a matter Lagrangian density, C^'^'>{4>) on each of the branes. 



2nl 



M 



d'x^\R--id^f-KlVi^) 
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J2 J d^xV^}^2 



«2c/W(0) + «2£W(^) 



(12) 



where M is the 5D manifold with metric gab, and the 4D hypersurfaces are the branes each having a projected 
metric h^^^. The constant Kg , which is effectively the 5D gravitational constant, has mass dimensions of 0{m~^) (m 
is a unit of mass), (p is dimensionless, the bulk potential V{(f)) has dimensions 0(rn^) and the brane potentials U'-'-\(j)) 
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have dimensions 0{m'^). The 2 [iiT^*^]]^ terms added to action on the brane are the appropriate Gibbons- Hawking 
terms 1 32 , required for a well-formulated variational problem. 
Varying this action with respect to the metric gives 

+ 2^ E d^xV^ I [Khab - Kabt + ^ (ui<jy)Kb + mhab - 2^) I Sg'^" , (13) 

where we have now dropped the i labels on the terms in the brane action for convenience; suffice to say, all the 
subsequent expressions apply on each of the branes. The [Khab — Kab]^ term is a combination of the boundary terms 
from the variation of the bulk Einstein-Hilbert term and that of the Gibbons-Hawking term (see ref. for more 
details as to its derivation). By demanding these expressions be true for all variations, one can read off the bulk 
Einstein equations and the Darmois-Israel junction conditions^ to be 

Gab = l:da^dbcl> - i (l{d<j>f + 4v{(t>)] gab , (14) 



2 ^-r "-r 2 \2 

[Khab - Kab]t = (^Tab " lu{cf>)hab^ , (15) 

where we have defined Tab, which we will interpret as the energy-momentum of 'ordinary' brane matter, by 
Varying (|l^) with respect to (j) gives 

where is the covariant 5D wave operator g'^^Va^b- Again, the i labels on brane terms have been dropped for 
convenience. Reading off the equation of motion in the bulk and the jump condition at each brane gives 

By using (^ and an appropriate generalization, one can write V'^^ in terms of the 4D covariant wave operator, V^, 
acting on (j) and normal derivatives of as follows 

= V^</) + Kn.d(t> + n.d{n.d4>) = K5 ^ ■ (19) 



This is expression is the result of the equivalent procedure used to derive (|1C|) for the Einstein tensor. However, as we 
shall discuss in the subsequent sections, it is impossible to write this equation in terms of local quantities only defined 
on the branes due to our lack of knowledge of n.d{n.d(f) on the brane. This is the main obstacle to the derivation of 
a closed form solution when there is a bulk scalar field in an equivalent way to the RS type models. 



C. Effective 4D equations 



Having derived the appropriate bulk Einstein equations and equations motion for (p along with their appropriate 
junction conditions, we are now in a position to derive effective 4D equations by applying the geometrical projection 



^These are almost universally referred to as the Israel junction conditions [Q, but this extrinsic curvature method was used 
much earlier by the French mathematician George Darmois |34|. We would like to thank Brandon Carter for pointing this out 
to us. Even earlier work on junctions in spacetimes was done by Lanczos |3q] and Sen M. 
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technology of the section [I A . The first step in this direction is to impose Z2 symmetry which relates quantities at 
each side of the brane and allows us to evaluate Kah and n.d(\) close to, but not on, the brane using the junction 
conditions, 



Substituting the expression for Kah^ and the bulk Einstein equation ( |14| ) into (|10|) gives the effective 4D Einstein 
equation. 

Gab = -A /lab + TkU{(t>)Tab + t|C/ T-Zlafc + Ks^^afa + ^ T" hab + 7;^ a4>V bcj) ~ TTT V^f /lab - Eab , (21 

12 lb oq) 61 \o(p J 6 24 

where, for convenience, we have written 

m = |c/(0)^ + ^V{4>) - |c/'(0)^ , (22) 

T^ab = --TTa'^Uc + TTTTTafc + -r'^'^Tcdhab — —T^hab ■ (23) 
4 Iz o z4 

A((/)) is an effective cosmological constant type term, nat is the quadratic correction to the Einstein equations first 



deduced, but interpreted somewhat differently, in ref. |11|, and Eab is the only non-local term, it being a projection 



of the bulk Weyl tensor onto the brane. In the next section we shall discuss how its effects can be understood simply 
in terms of quantities which are local to the brane in the case where the metric on the brane is spatially isotropic and 
homogeneous. 

By taking the trace of (^ij) we can obtain an expression for the 4D Ricci scalar, 

B = 4A(« - - 4 p - 3.-..,) - |(^^)% \m' . (24) 

which is independent of the contribution from the bulk Weyl tensor Eab and will be useful when considering the 
expansion rate. 

Similarly, substituting (|2^) into the Codacci equation (^) on the brane gives us the following expression for the 
divergence of the brane energy-momentum tensor 

b 1 S /L 

vVah--2^a(/'^, (25) 

which is an effective energy momentum conservation equation for the matter on the brane. Interestingly, the brane 
energy-momentum is only conserved either if C is independent of (p, or if cj) is constant on the brane. 
One can also substitute (^o|) into the the equation of motion for (p on the brane which gives, 

V'^ + n.din.dc/,) = KlV'icP) - ^{r - 2U{cf,)) (u'{cf>) + . (26) 

We have written this in terms of a d'Alembertian wave operator on the brane, local source terms, and the second 
derivative of 4> in the direction of normal to the brane. This term, which is analogous to Eab in the case of the effective 
Einstein tensor, is the only term which is dependent on the bulk. Unfortunately, in contrast to Eab, it is not possible 
to understand it simply in terms of quantities local to the brane and hence in our subsequent discussion we will be 
forced to make some assumptions as to the behaviour of on the brane. In previous discussion we have pointed 
out that RS showed that the theory of gravity on the brane was a tensor theory, and we see that in the effective 
4D Einstein equations are only slightly modified by the bulk in the weak-field limit, but there is clearly no way of 
thinking of the bulk scalar field in the same way; it being a truly 5D quantity, and must be treated as such. 

Now, motivated by the isotropy of the CMB and homogeneity of the observed galaxy distribution, we choose the 
metric on the brane to be an embedding of the 4D FRW universe, that is, 

habdx°'dx'' = -dt^ + a{tf-f^jdx'dx^ , (27) 

where ^ij is the metric of the three-spaces of constant curvature. In this case, standard results give 
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iCoo = + 4 ' (28) 
R^&{h + 2H'^ + ^ , (29) 

where dots denote difFerentiation with respect to cosmic time, a is the scale factor and H = a /a is the Hubble 
parameter. Therefore, in order to derive equations for iJ, and hence an analogy to the Friedmann equation, we have 
two alternatives: either compute Eqo on the brane and hence Gqo as advocated in ref. |0, or integrate the equation 
for H in terms of R which implies the existence of an integration constant. This exactly the equivalent approach to 
those taken in refs. ||ll|,|l5|. In the subsequent sections on applications we shall perform these calculations explicitly 
for a number of cases, to show that they are equivalent. For the rest of this section we shall attempt to compare our 
results with particular 4D theories of gravity using the later of these two approaches and in the next section we will 
discuss the precise dynamics of Eat on the brane. 

If we now assume that C is independent of (j) and that the quadratic corrections are small, then 

i? = 4A(0)-^C/(0)r + i(W)2. (30) 
The equivalent expression for Einstein gravity with a cosmological constant is 

Gab = SnGNTab - Ahab , i? - 4 A - SttGnt . (31) 

Therefore, we see the directly the interpretation of A as a cosmological constant and if one defines the an effective 
gravitational constant given by 

^C/((/.) = 8^G^(</)), (32) 

then one we will get the same expansion rate as Einstein gravity in this limit. As we have already discussed, from the 
point of view of the brane we have no control over the dynamics of the scalar field, but the fact that it is related to 
the effective gravitational constant means that it must not have changed much since the time of BBN in any realistic 
model of the universe. Note that in making this definition we have justified the assumption of ignoring the quadratic 
terms at late since they will be suppressed by a factor of t/U. 

Clearly, by making these identifications we can make some interesting general statements. First, we should note 
that the expression for the Ricci scalar is independent of Eab and hence, if <j) is constant on the brane, it is completely 
insensitive to anything which is not on the brane. This is not to say that the gravitational dynamics are in general 
independent of the bulk, just that in this specialized case the only effect on the effective Friedmann equation is via an 
integration constant. Secondly, we see that U{(j)) must be a strictly positive number, at least after the epoch of BBN, 
in order to recover the desired expansion rate. We will discuss this in detail in the context of the Goldberger-Wise 
mechanism for stabilizing the extra dimension; suffice to say if the effective tension on the brane is negative then 
the corresponding expansion rate will also be negative. This is effectively rules out the RSI model for solving the 
hierarchy problem, even when there exists a scalar field to stabilize the extra dimension. Finally, in such models there 
is a simple and at first sights uncomplicated way to ensure a zero effective cosmological constant on the brane, by 
setting 

t/'(0)'-|^(</')' = (33) 

an equation which can be solved analytically in the simple case where V{(t>) = and V{4>) is a constant. This self- 
tuning mechanism, first suggested in refs. |3^j3^, has been the subject of a great deal of debate in the context of 
supergravity models where U and V are derived from the same super-potential (see also [ ^9po| ). Unfortunately, it 
appears [ p4| that no such supersymmetric model exists without the a naked singularity somewhere in the bulk. 

We have already noted that it was deduced in ref. that there is effective BD gravity on the negative tension 
brane in the RSI model due to a localized mass distribution on the positive tension brane, and this may at first sight 
appear to in conflict with our results. However, here we are only considering the effect on expansion in the spatially 
homogeneous and isotropic model, and hence their assumptions do not apply. In fact, a localized mass distribution 
off the brane cannot contribute the overall expansion rate of the brane since its only gravitational effect must be via 
Eab which is traceless and hence does not contribute to Ricci scalar. 

Even though the analysis of ref. docs not apply here, it is sensible to compare our model to alternative theories 
of gravity such as BD, rather than just Einstein gravity since it can help us to constrain further the parameters of 
our model. The field equations for BD gravity in 4D without a cosmological constant are M 
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-Tab 



if if- 

and the BD field, Lp, has equation of motion 



^a'^bf - Kb^'^f 



—2 

V = 



87r 



3 + 2tj 



Here, we have kept the bars and hab as the 4D metric to avoid ambiguity. Taking the trace of 
(^) gives an expression for the Ricci scalar 



and hence we make the following identification 



2w 



3 + 2cj 



Stt 



2u) 



3 + 2w 



(34) 
(35) 

|) and substituting 
(36) 

(37) 



Thus, if U{4)) < and ip > 0, the Brans-Dicke parameter, cu, must be in the range —3/2 < ui < 0. Note that the 
same bounds on w were found in pO| ] when considering localized mass distributions on the negative tension brane. 
Such values of to can be ruled out by experiment reiterating that the RSI model, and associated models with a 
stabilized scalar field, are ruled out by experiment. 



D. Conformal Transformations 



There has been considerable debate in the literature of scalar-tensor gravity theories (such as Brans-Dicke theory) 
as to which frame represents the physical one. This issue is also pertinent in the case of dilatonic brane-worlds due 
to the existence of a scalar field. Scalar-tensor theories can be related by conformal transformations of the metric, 
gab — > dab = ^{4')^9ab and it is well understood that changing to a different frame in a cosmological solution can 
change the expansion rate (in particular, it could change the results of the following section. We have adopted the 
view of that the frame in the formulation of the action ( |l2|) is the physically significant one in the sense that 
it is with respect to this metric that free particles will follow geodesies. As this viewpoint is not universal, even 
in the 4D case (see for a summary of different opinions) it is interesting to study how our effective Einstein 
equations behave under conformal transformations of the metric. In particular, we can demonstrate that, unlike with 



the Jordan frame of Brans-Dicke theory |48 49[|, the frame chosen by the induced metric is the one in which the 
energy-momentum of the scalar field does not contain second derivatives of (j). This means that any other frame would 
be philosophically problematic in that second derivative terms necessarily give an energy-momentum which is not 
positive-definite. (Again, there are differences of opinion: not everyone considers violations of positive-definiteness in 
the gravitational scalar sector to be problematic ]4^!) 

If the 4D induced metric transforms as hab — > hab = ^{4')'^ hab then the covariant derivatives which preserve these 
metrics are related by, for example, 

Vaffc = VaVb + \habh'"^dd(t) - 5l^(t)d^(j}j Vc , (38) 

and the Einstein tensor transforms to (see p^) 

Gab = Gab - 2^VaV6(/) + ^ (20'' - 0") V^^Vfe./) + ^ {ft'V ,V + fl'V ,0 d^) habh'"" ■ (39) 

Note that ria = and n" = O^^n'' and hence Eab must be conformally invariant as both C^bcd and n"nh are. 
Likewise, and n'^b are conformally invariant. Substituting for Gab from ( pl| ) and expressing everything in terms 
of derivatives relating to hab, one can deduce the effective Einstein equation in the new frame, which is 

Gab = -^hab + jl^Tab + ^^-ab - Eab ^ 2-V.V,0 + (^3 + 2^ " 2^ j V.0V,0 

^^~^?^~l]m'}hab. (40) 
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Note that we can substitute for V^^ using the conformal transform of (J^^ and that some of the terms can be set to 
zero by suitable choice of f2. However, the VaVb0 term will always be present unless is constant, which, of course, 
corresponds merely to a rescaling of units. Thus the original frame is the only one in which the terms corresponding 
to the effective 4D energy-momentum of the scalar field can be positive definite. Also note that, in all frames, the 
coefficient of Tab has the same sign as the vacuum energy, as in Brans-Dicke theory. Frames other than the one we have 
considered have an additional term proportional to rhab', the coefficient of this term is constrained by gravitational 
tests such as those carried out on solar system measurements. 

In order to consider the expansion rate in the new frame can also calculate the equation of motion for the transformed 
Ricci scalar, which is 

R = ^A + ^r + ^n+i^4-[^\+-^[-[^\\ +6l^V^0+|6:- + 14l--^)(V0)^ (41) 





As mentioned above, it is possible to substitute for V2(/) using (^6|), which gives 

= 2^(W)' - n.din.dcf.) + (,/>) " ^(r - 2C/(0)) (^C/'(</)) + ^ j j . (42) 

Thus we can remove the second brane derivatives in </> at the expense of introducing the second derivative normal to 
the brane, n.d{n.d(f). This gives 

i? = 4A + -I— + -2—— T 





4 + 2^-l)(v*)^-6^^„,a(..a«. (43) 

where we have defined 

1=^ + ^{&4V' + 4UU'). (44) 

It can be seen that the coefficient of r in this frame is different from the coefficient of t in the original frame. In 
particular, is is possible to choose ^{(f) so that the expansion rate is positive even when the brane has negative 
tension. However, the effect of the n.d{n.d(j}) term is impossible to quantify. From the expression for the Einstein 
tensor, it is possible to see that, although one can change the sign of the dependence of the Hubble parameter on the 
matter content by changing frame, one cannot change the sign of the effective Newton constant. 

Our original choice of frame for the 5D action was the 5D Einstein frame and one might want to argue the case for 



other 5D frames. The result should be recoverable from our conformally transformed equation (40). This has been 
studied in |50[| starting from a 5D action in a general frame: 

S=-L f d^^V^{F{^)R + ---}~^Y. f d^^V^\2Fi<t>)\K^^A'^ + ■■] . (45) 

We refer the reader to for a discussion of the results and merely quote that the effective Einstein equation derived 
from such an action has the form 



Gab = SnGNTab + S-KGAThab H , (46) 

where Gn is in agreement with our conformally transformed equation ( ^ ) above as expected. In particular, it should 
be noted that Gat always has the same sign as the vacuum energy of the brane. 

In summary, although there are differences in opinion in the literature of scalar-tensor gravity theories, the frame we 
have used is in good agreement with the philosophy of both |^] and , in the first case because we have interpreted 
the metric of the original frame as that governing the geodesies and in the second because ours is the only frame in 
which the scalar field has positive definite energy-momentum. 
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E. Understanding the effects of the bulk 



Our 4D effective equations for the gravitational sector and the scalar field (|2^) both contain terms which specify 
how effects from the bulk can manifest themselves on the brane. In the case of the scalar field it is impossible to 
derive a simple expression in terms of just local quantities forcing us to make sensible assumptions as to the behaviour 
of (p on the brane. However, it is possible to do so for the in the gravitational sector. 

One can deduce an equation for Eab on the brane by realizing that the 4D Einstein tensor must satisfy the 4D 
Bianchi identity V Gab = 0. From one can deduce that 

oq) lb 0(p 

* i (^'<« + 1) (i) + 5^'*^'* - ii^'<^*'' ■ <«' 

The tensor Eab is traceless, tangential to the brane and satisfies the above equation. Therefore, from the point of 
view of the brane it includes 5 independent quantities, 3 of which are related to spatial isotropy of the brane via E^i 
and the other 2 are tensor degrees of freedom due to gravitational waves in the bulk. 

If we assume that the brane and the scalar field on the brane are spatially homogeneous and there are no bulk 
gravitational waves then E^i = since Goi = 0, and Eij — f{t)^ij for some function f{t). Therefore, Eab is entirely 
specified by the 6 = component of (|47| ) , 

i?oo + 4^i?oo = |a'((/.) + BirG'MVm + 47rGAr(</.)^ - ll^"^'^)^ " i*^} " ^^^"^^o (48) 



16 [ Scf)] dt _ 

If C is independent of (j) and we consider a perfect fluid, we can write T°'b ~ diag(— and hence we find that 
V ttqo = 0. Thus 

Eoo + 4^^oo = ^ |a'(</.) + 87rG^(</.)p - ^4 , (50) 



and so, if (p is constant. 



l|(a^i^oo)^0. (51) 



which integrates to give i?oo = Ma~^ where the integration constant, M, has dimensions 0{w?). This is the exact 
same integration constant as discussed in terms of integrating the equation for the Ricci scalar in the previous section. 
Clearly, if (p is slowly varying then modifications to this result can be derived in a perturbative expansion. 

We have shown, therefore, that the effect of the bulk on a spatially isotropic and homogeneous brane is to contribute 
an effective radiation type term to the Friedmann equation ||l^-|l^,^,|^ . The physical interpretation of the integration 
constant is that it quantifies the mass outside the brane, which in the case of RS2 is zero and in RSI is just that of 
the other brane. This is almost equivalent to Gauss' law in electrostatics and is the analogue of Birkhoff's theorem 
for spherically symmetric space-times (see [Q for a discussion of Birkhoff's theorem applied to domain walls in 
adS space). In fact this result is much deeper. In ref. it was shown that the only global spacetime compatible 
with spatial isotropy on the brane is Schwarzschild-anti-de-Sitter and this integration constant is the equivalent of 
the Schwarzschild mass term in usual treatment of vacuum spacetimes in 4D. Therefore, the precise magnitude this 
contribution to the expansion rate not only depends on the mass of other branes, but that of any black holes nucleated 
by thermal effects in the bulk p6| , ^ . Since it behaves like radiation on the brane, the magnitude of such a contribution 
to the expansion rate can be constrained using BBN. 

A corollary of this, reiterating our earlier point, is that when <f> — Q, there is not possibility of changing the sign of 
the effective gravitational constant due to effects from the bulk. Hence when p is stabilized a negative tension brane 
will have a negative expansion rate. We should note that this does not preclude such effects when the scalar field 
is varying, nor does it prevent effects being transmitted between branes by perturbations from spatial isotropy and 
homogeneity, for example, the localized mass distribution of ref. ]To| . Such effects would, however, be at perturbative 
order in the cosmological context. 
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III. APPLICATIONS 



A. Randall-Sundrum type models 



The basic RS scenarios are a very special case of the general class of models we have discussed where is a constant 
at all times. The potentials will then take constant values, related to the CCs by 



(52) 



where A is the bulk CC and are the brane CCs with i = 1, 2 and Ai = — A2 in RSI. Substituting these potentials 
into ( ^ ) gives the result obtained in ref. |l^], that 



Gab — [ 1 Y' j ^<^^ ^ 6~^°^ ~^ '^h''^ab — Eab , 



(53) 



and the standard energy conservation equation V Tab = 0. Note that this result holds for both the RS scenarios: 
there will be an equation for each brane in the first scenario, but that the i labels have been dropped for convenience. 
From this we can deduce the the only way in which information about the bulk can be communicated to the brane is 
by the Eab term. The first term, which is the effective 4D cosmological constant, is usually set to zero by tuning A, 
but we shall not assume this to be so here. Note, once again, that recovering a positive expansion rate requires that 
A > 0, as found in ref. p7[|, and on our brane, we identify 



G 



N 



4x 

487r 



(54) 



We can recover the Fricdmann equations for these solutions using our formalism in two ways as described in the 
previous section. Using (^8|) and assuming that the matter on the brane to be a comoving fluid with energy-momentum 
tensor given by T°f, = diag(— p,p,p,p), one can deduce that 



4A. 
~6~ 



SvrG 



N 



^5 2 1 ^ 
— p - -Eqq 

36 3 a'^ 



(55) 



This is the now well known Friedmann equation for a RS type model and includes a term analogous to that found in 
4D, a quadratic order correction and the often ignored unknown Eqq. When the brane metric is spatially isotropic 
and homogeneous, we can determine Eqq from ( |5l| ) which give 



4A 



36 



SttG 



N 



Kt , M k 



(56) 



The other way to recover this solution is to use (|29|), which gives 

AttGn 



4A 



18 



3pp) 



(57) 



in combination with the conservation of energy-momentum on the brane, p + SH{p + p) = 0, which can be used to 
remove p in favour of p and p. If the dependent variable is changed from t to a, the differential equation for 
derived in ref. [ nSl is recovered and the solution (56) is once again found. In fact, this analysis is slightly more general 
than that of refr||l5) as it only requires homogeneity and isotropy on the brane, not necessarily in the bulk. 

One clear problem of these simple RS type models is the fine-tuning of the CCs when one tunes A to make the 
effective 4D cosmological constant zero to get a Ricci-flat brane in the absence of matter. It has been suggested that 
this unwanted aspect of the models can be removed by using a scalar field to stabilize the extra dimension in an RSI 
type model, and our a ppr oach has been set up to deal with this kind of scenario. Such a mechanism was proposed by 
Goldberger and Wise with 



vi<p) 



2A. 



fJ-i 



2A, 



(58) 
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where the fj, and /ii have dimensions 0{m), while the Vi are dimensionles^. If the field is varying, the equation for 
Eqo can no longer be integrated exactly, nor can the analogue of (^7|), both due to the rather complicated dependence 
on (f>. The expansion rate will in general be very different and hence the scalar field must stabilize before BBN. Not 
surprisingly, however, if the field is constant then the whole situation just returns to that of a simple RS model and 
the integration can once again be done exactly. If the field is slowly varying then, at zeroth order, the solution will 
be that for </> constant with small pcrturbative corrections. The model stabilizes when the takes the values ±i>i on 
the branes. Thus, after stabilization, the coefficient of p in the expression for is the same as (^6|). Hence, if Ai is 
negative then negative expansion is expected at late times as in the simple RS model. 

What we have deduced in a rather roundabout way — adding the scalar field and then it making constant after 
stabilization — is that the scalar field can have no effect on the sign of the expansion rate once it is stabilized; the 
reason being that, as (^) comes from the trace of (^l]) and Eab is trace-free, there is no contribution from the bulk 
Weyl tensor in (|57|). Thus the only contribution to that can possibly arise from the bulk Weyl tensor is that 
term arising from the integration constant when solving (p7[). This is the argument discussed in the previous section 
on the general case applied to the specific case under consideration here. In particular, we note that recovering the 
standard Hubble expansion when p is small req uires that A be positive^. The same result was discovered in jsill , where 
an explicit 5D metric was used. From section II D we see that it is possible to recover a positive expansion rate for 
cosmological solutions in the case of a negative tension brane, although this is unsatisfactory as it does not give the 
usual Einstein equation and leads to an indefinite contribution to the energy-momentum from the scalar field. 

This result is almost the total opposite of that obtained in ||l8|-^. In these papers, the standard Friedmann 
equations are recovered from the RSI scenario by averaging over the bulk to obtain an effective 4D theory in a similar 
way to the dimensional reduction performed in Kaluza-Klein theories. In simple terms, these authors make a specific 
coordinate choice and integrate over one of the coordinates. To compare this to our covariant formalism of projection 
onto the brane, we now attempt to formulate such an averaging for a general manifold with co-dimension one branes 
without reference to any particular coordinate system. 

The covariant analogue of integrating over the bulk coordinate at each point on the brane is to integrate along a 
congruence of curves. Note that to integrate tensors we must use a pull-back to bring them into the tangent space at 
the point on the brane. This averaging procedure will depend on the congruence chosen and, in the case where there 
is more than one brane, on the brane chosen. For any brane, the natural choice for the congruence is clearly given 
by the geodesies tangent to the normal on the brane. This gives rise to the same foliation discussed in section ||, 
which is well-defined near the brane but may not defined throughout the bulk. The most serious objection to this in 
models with two branes is that the geodesies may not be normal to the second brane and so a different result would be 
obtained working from the other brane. Note that the method of projecting onto the brane only requires the foliation 
to be well-defined in a neighbourhood of the brane. 

This brane-world averaging should be contrasted with the usual Kaluza-Klein dimensional reduction where there 
are no branes and all matter experiences all of the dimensions. (Of course, our projection method is meaningless 
in theories such as these where there is no preferred hypersurface.) In such theories, everything is considered to be 
relatively homogeneous across the extra dimensions, i.e., there are no preferred points. In contrast, with brane worlds 
the matter has support only on the branes and, in addition, the background is usually highly warped. The effect of 
integrating over the bulk with the exponential weight factor is to give 4D quantities which are largely determined 
by the matter on or close to the shadow brane. It is questionable whether it is sensible to interpret the average as 
physically meaningful, especially as our experiences relate to a very atypical point in the spacetime. 

As the Goldberger-Wise method for stabilising the radius involv es in troducing a scalar field, it is possible to re- 
interpret the results by a conformal rescaling, as discussed in section II D . Some authors, such as ||2^, use a conformal 
scaling chosen by integrating over the bulk. For the cosmological solutions they consider, this gives a positive expansion 
rate, even for negative tension branes: an effect which can be seen by considering the term proportional to r in (|4^). 
However, as we have already observed from (^) this does not recover a positive Newton constant in the full Einstein 
equations. 

The difference in the results obtained by averaging and projecting raises the philosophical issue of which should be 
interpreted as the correct 4D effective theory. This difference does not seems to have been discussed in the literature, 



*We should note that our definitions are such that the dimensions and the particular forms differ from those used in ref. [ p2[ . 
However, the physics are exactly the same. 

^One could recover an expanding universe with negative A during the radiation era by choosing M sufficiently large to dominate 
the ordinary radiation, but one would never get matter domination, since as soon as the matter terms begin to dominate the 
negative expansion terms would kick in and contraction would ensue. The process of expansion and contraction would take 
place indefinitely in such a model. 
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where some authors, e.g., p^-|2l]], have adopted the averaging method whereas others, e.g., | |Tl| , p^ , ^ , have obtained 
different results by hands-on versions of our method where a specific form of the 5D metric is used. In particular, 
concludes, as we do, that the brane tension must be positive to get the correct cosmoiogicai expansion even in 
presence of a stabilized scalar field, whereas |^ conclude the opposite by an averaging analysis. This discrepancy for 
cosmoiogicai models could be resolved by reinterpreting the metric via a conformal scaling, but this is problematic 
for reasons we have already discussed. One technical advantage of the projection procedure is that it can be defined 
on a larger class of manifolds. On a more philosophical level, the defining property of brane-world models is that we 
are confined to the brane and so the only matter we observe via non-gravitational interactions is that on our brane 
and the curvature quantities we measure are those of the metric restricted to the 4D submanifold, so the projection 
method seems a more natural interpretation. 



B. Hofava-Witten model 



The Hofava-Witten model is a realistic scenario for a brane-world based on the a compactification of IID description 
of Eg X Eg heterotic string theory, often called heterotic M-theory The general model contains many fields, most of 
which can be set to zero self-consistently. However, one, a scalar field, cannot since it corresponds to the deformation 
properties of the compactified Calabi-Yau manifold which is used to reduce the model from IID to 5D. In this section 



we shall explore the implications of similar models with a scalar field. As a result of our analysis in section IID, we 
will work in the 4D conformal frame selected by the induced metric. 

Within the general framework in which we have been working, we shall consider two branes and the potentials 

- «5-^°/'Ce-2^* , C/,(0) = K-'^'Ce-^'^ , (59) 

where the powers of K5 have been chosen to make the C and Ci dimensionless. We expect these potentials to 
correspond to 'slow-roll' behaviour of (f). Spherically symmetric cosmoiogicai solutions for dilaton domain walls with 
these potentials have been found in ref. using Birkhoff's theorem. The effective Einstein equation now has 

-4/3 4/3 
A(<^) = ^e-'^^ (24C + (2 - 3(3^)Cf) , 8^Gjv (0) = '^Qe~^'^ , (60) 

which leads to a variable cosmoiogicai constant term, similar to quintessence, and also a variable effective gravitational 
constant. From ref. jsj we see that, for the simplest HW theory, /3 = 1 and 24C — Cf; thus A = in this model. We 
should note, however, that the mechanism by which the a flat brane (A = 0) is achieved is subtly different from the 
RS type models since the potential in the bulk is positive; the cancellation being due to spatial gradients of the scalar 
field. We also have Ci — — C2 so, from the previous analysis, we must live on the brane which has a positive value of 
Ci to recover Einstein gravity. 

We study the cosmology of this model using the same methods we employed to the RS type models. The equation 
for H (^9|) gives us the following differential equation: 

4/3 4 11. 

H + 2H' = \c.e-^\p - 3p) - + 3pp) - ^^'^^ - ^ , (61) 

where we have assumed, for simplicity, that C is tuned in terms of Ci and /3 so that the cosmoiogicai term vanishes 
identically. The equation for E'en (Esl) becomes 



4^i?oo = ^0 (3?i - 4l^pC,e-P^^ , (62) 



and we have the standard energy-momentum conservation condition p -\- •iH{p + p) = 0. As we have already noted 
the effective gravitational constant varies, but in the way which is different from that usually assumed in 4D. The 
standard way to allow Gat to vary in 4D GR is to make V (GNTab) = as a consequence of the Bianchi identities, 
but this is no longer true in this case and Gat can vary without a modification to energy-momentum conservation. 
Hence this HW type model provides a novel mechanism for the variation of Gat 

At late times, we will require (p to be almost constant. The solution of zeroth order in (j> is effectively the RS 
solution, but in a specialized case we can attempt recover the next-to-leading order correction in (jj. If we ignore the 
(jx^ term, the equation ( |6l| ) still cannot be integrated in general due to the (f) dependence in the first term. However, if 
the matter satisfies the equation of state p = 3p, that is, if the universe is radiation dominated, then the first term will 
be identically zero. Integrating the energy conservation equation, gives the usual solution for radiation, p = p^a^^. 
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where Pais the density of radiation today, if we assume a(to) = 1 and to is the time of the present day. Thus we can 
write (|6l| ) as 



dH^ .rr2 f^tpl 2fc 



which can be solved to give 



where C is a constant of integration. This integration constant corresponds to the combination of the bulk effect and 
the effect of the radiation localized to the brane. 

We can also attempt to derive this solution by computing Eqq. With p — pf^a~^, ( |6^ ) can be solved to first order 
in derivatives of 0; 

- E„^l,{^^^C,e-^*-^y (65) 

where we have chosen — 3/C as the constant of integration for consistency with (|6^. Equation ( p8| ) gives us 

2 ^ >^^c,e-P^ + ^ _ 4 - -^^00 . (66) 
36a4 36a8 3 ^ ' 

We note that, when we substitute the solution for i?oo: there is a cancellation between the terms dependent on (f) 



and the result is consistent with ( 64 ) . Thus we see that there are no first-order corrections in derivatives of to the 
Friedmann equation in the radiation dominated era. There is, however, a 4> dependent correction to E'en which cancels 
with the term proportional to p in (|6^). We conclude, therefore, that the RS type solution is valid if one just ignores 
term which are second-order in derivatives. 

This solution is valid during the radiation era and when is rolling slowly enough for any second-order terms in 
derivatives of (j) to be neglected. By the time the universe has reached the matter domination, one might hope that 
the evolution of (j) is now sufhciently slow that (p will be almost constant, so that all derivatives of (p may be neglected 
and hence the RS solution can recovered in the matter era as well. Thus one might imagine that the cosmological 
solutions during matter-domination and during radiation-domination, except possibly at very early times when is 
not changing slowly, are those of the RS model. Given the very different nature of the HW and RS models, this is 
somewhat remarkable. 

C. Cosmology with matter coupled to (f) 

Our setup allows us to study the case where the matter on the brane is coupled to the scalar field 0, i.e., when 



SC/6<j) ^ 0. We now study such a scenario in a cosmological context. Since (24) has SC/6<j) terms, we need an explicit 
form of the Lagrangian. Since we want a cosmological solution, we want the energy-momentum tensor on the brane 
to be Tab — F{4>) {phab + {p + p)UaUb): F{(j)) is an arbitrary function, p is the energy density, p is the pressure and U 
is the velocity of the flow. This suggests that we try F{(j)) times the Lagrangian for a perfect fluid^ so we consider 
the variation where the matter part of the action is given by 

^matter = -2 j ^d^xF{cj)) [p{e, s) - ^ (/I'^'f^.f^b + E^)] . (67) 

Here, s is the entropy, e is the enthalpy, n is a Lagrange multiplier and f2 is a 4-vector field, which can be written 
as rta = VoX + ctV a(3 + OV aS for some scalar fields Xi ol, /3 and 6: respectively the three Clebsch potentials and the 
thermasy. Varying ( pT]) with respect to n gives the constraint h'^'Tla^b -|- = so we define the velocity U — Sl/e. 
Varying with respect to the other parameters e, Xi ^ and s gives 



^See [p2|l for details of how this is constructed. 
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^=n, Va[nF{(b)U'']=0, aP ^ aU ^ aS ^ Q , ^+nU''Va9^0. (68) 

oe OS 

Varying with respect to the metric gives, as required, 

Tab = F{^) [phab + neUaUb) = [phab + [p + p)UaUb) ■ (69) 

Thus, on-shell where If^Ua = —1, the Lagrangian and its variation with respect to cj) 

m = -2F(0)p(e, s) , ^ = -2F'(0)p(£, s) . (70) 

From (|l|), we see that the effect of couphng the matter can be thought of as changing the relation between p and p. 
CoupUngs of scalar fields to photons and baryons are tightly constrained; coupling to dark matter is less constrained 
(see [p3| and refs. therein for more details). A s befo re we have chosen to work in the frame determined by the metric 



induced on the brane. As discussed in section II D, there is an issue of which frame to choose. An alternative to the 
above approach would be to transform the action into the Jordan frame before varying the action and then transform 
the equations of motion back to the Einstein frame. The variation of the action in arbitrary frames is given in [50| ]. 

Now that we have the expression for the Lagrangian that we need, let us turn our attention to the cosmology of 
such models. If we assume that V{(t)) has been tuned so that A = and that 4> is spatially homogeneous on the brane 



(71) 



(i.e., (f> on the brane depends only on t) then ( |24| ) and ( ^91 ) give us 

H + 2H'^ ^Ummp 3p) + ^U'{^)F'{^)p - ^Fic^fip' + 3pp) ^F'(0) V - ^'P' ^ ^ 



and (|25|) gives us 



p + 3H{p + p) = -^^{p + p). (72) 



into (|71|) gives 



For an equation of state p = wp, ( ^2[ ) has solution p = po (a'^F) Por dust {p = 0), substituting this solution 

iJ + 2H-^|w)^-fM-l^^-i, (73, 



which is identical to the non-coupled case. This is not surprising as, the energy-momentum tensor is r^f, = 
diag(poa~'^i 0, 0, 0) so we have the same matter content as for non-coupled dust. It is also a promising start for 
this model as the standard Big-Bang model has been very successful as a description of the matter-dominated uni- 
verse and so we want our model to look like the standard cosmology during the matter era. 

Now consider radiation matter in the specific case where </> is constant. For convenience, we will assume that p^ 
and (j) terms can be neglected. Then (^l|) becomes 

^{a*H^)^A-~2kad, (74) 
at a 

where A is a collection of constants. We get the solution 

For the non-coupled radiation matter, A in ( [Tq ) is always zero but the coupling allows a non-zero, positive value. As 
a grows, this new term will dominate the a~'^erm. Also, as a — > 0, the first term will dominate and will be negative. 
Thus there will be a value of a, dependent on the values of po, C and k, for which = 0. 

The term proportional to a~'^ contains a contribution from the bulk and one from radiation matter on the brane. 
To calculate these individual contributions, we can use ( p8| ) which, ignoring quadratic effects, gives 

= + u'F'F-^/') Ieoo + ^ . (76) 



We can calculate Eqo from (49) which reduces to 
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^(«i^oo)= 3^,/3 ^(loga), (77) 

where quadratic matter terms have again been neglected. This generates the a^^ log a term and the other part of the 
a^^ term in ([75[). 

We have a Friedmann equation which is like that of the standard cosmology but with the extra a~^loga term. 
Motivated by this, we now study the cosmology generated by a Friedmann equation with this new term in addition 
to the the usual radiation, dust and curvature terms. Thus we consider 

^ Aloga B C k , , 

where A, B and C are positive constants. We will consider this very much as a dynamical system, although we will 
make some physically motivated assumptions about the relative magnitudes of the constants. 

To determine whether this models suffers from a flatness problem, we define the curvature density in the usual way 
by Vlk = —k/{a'^H'^), whence we derive the relation, 

nk^-^{H + H'^nk- (79) 
This tells us that fife = is a stable fixed point if ij + H'^ > 0. From (Fsh, we get 



So there will be value of a, dependent on A, B and C, below which ilk = will be an attractor. This stability 
condition can be written, 

A-2B-2Aloga-3Ca>0. (81) 

The r.h.s. of ( [78| ) can be negative for some values of a. Clearly we must have > 0, so we also require that a only 
takes values for which the r.h.s. of ( |7^ ) is non-negative. We will neglect the last term both because its relative effect 
is only significant for large a and because we hope to ensure that the universe is driven towards flatness. Thus we 
must also satisfy the condition 

A\oga + B + Ca>0. (82) 

If we can satisfy both (|8l|) and (|8^) we will get a period of accelerated expansion which drive the universe towards 
flatness. Adding twice the second inequality to the first gives us A > Ca. Substituting a — A/C into the second 
inequality gives 

,„g(5)<l + f^ (83) 

If this is not satisfied, then ( ^2| ) cannot hold for any value of a < A/C. Thus, the coefficients A, B and C must 
satisfy ( ^3| ) for us to get a period in the history of the universe where £7^ = is an attractor. Note that this is a 
necessary but not sufficient condition. We now consider whether we can achieve a period of accelerated expansion for 



A \ I A log a 3C 



some reasonable values of the parameters A, B and C. Define B — B/A and C — C/A. Then conditions (81) and 
( p^ ) become 

/i(a) = 1 - 2B- 21oga-3Ca > 0, (84) 

/2(a)=loga + B + Ca>0. (85) 

Now /i decreases monotonically from /i(0) = oo whereas /2 increases monotonically from /2(0) — —oo. Thus Ela > 
such that fi{a) = /2(a). If /i(a) = /2(a) > then there is an interval I 3 a such that, Va G /, /i(a) > and 
/2(a) > 0. 

The ratio C/B = C/B is well known and, if we set a = 1 at present times, it is 10'* to good approximation p4[ . 
By plotting /i and /2 against a for different values of B we can determine that, for B > 11, there will be a period of 
accelerated expansion. Choosing larger values of B makes this period occur for smaller values of a. In the standard 
cosmology, BBN occurs when a ~ 10^*^. So we expect the period of accelerated expansion to finish before this. 
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Choosing B = 30 gives a period of acceleration 10~^^ ^ a ^ 1.5 x 10~^^, which would be acceptable. In practice, B 
is probably much larger than this as no effect like the a"'' log a term has been observed. 

As discussed above, there is a finite value of a for which H = 0. Clearly, it is important that this is an unstable 
fixed point so that the universe expands from this initial value. Hence we now analyse the stability of the fixed point. 
Changing (^8|) to conformal time, we find 

-^l = A\oga + B + Ca-ka^ . (86) 
or J 

Scaling t ^ I3t and a ^ ja gives 



f}' \dr) 



2 



^ A\oga+ {B + Alog"f) + Cja- kYa\ (87) 



so we can choose (3 and 7 to give 



2 



a'2 = ( ^ ) ^loga + Aa-Ba\ (88) 

where A and 13 are redefined constants with A being positive and B having the same sign as k. Primes denote 
differentiation with respect to the rescaled conformal time. If fc < then a' = has exactly one root, whereas, if 
fc > 0, the number of roots will be two, one or zero depending on the values of A and B. The ratio of A to i? is 

I = je^/^ . (89) 
B k 

Since we expect B ^ A for a reasonable cosmology, this ratio will be large, so there will be two roots. The case k > 
with fewer than two roots docs not allow a'^ > 0, so does not produce any sort of cosmology. 

To analyse the stability of this fixed point (the one with the smaller value of a in the case where there are two) we 
calculate a" by differentiating (|8^. We find 

2a" a^^ + A-2Ba. (90) 

If fc < 0, a" > Va. If fc > then 3^ such that a > ^ ^ a" < and a < ^ ^ a" > 0. However, A/B > 1 so 
^ 3> 1 and so we have a" > at the fixed point with the smaller value of a. Thus, in all cases, we have a cosmology 
where the universe expands from a finite value of a. In the case fc < 0, the universe will expand indefinitely whereas, 
if fc > 0, the universe will expand to the other fixed point. To determine the stability of the other fixed point in the 
latter case, note that 

a" = - (- + -a- Ba'^] <loga + A- Ba"^ ff a > e^/^ (91) 
a \2 2 I 

Thus a" < at the second fixed point if a > e^/^, which is likely as A^ B, so the universe will recoUapse. 

The upshot of of all this analysis is that we have models where the universe expands from a finite value of a 
with an early period of accelerated expansion (i.e. where f2fc = is an attractor). These cosmological scenarios are 
interesting as they avoid an initial singularity and alleviate the flatness problem of the standard cosmology. In the 
above discussion, we have ignored the terms quadratic in the matter density. Thus the accelerated expansion must 
occur at low enough energy scales for the quadratic effects to be negligible but at high enough scales not to affect 
BBN. 



IV. CONCLUSIONS 



Using the Gauss-Codacci formalism we have derived effective 4D Einstein equations for a brane-world scenario 
with a bulk scalar field. By studying the trace of these equations we have been able to circumvent our ignorance of 
the contribution from the bulk in the gravitational sector, but not for the scalar field. This allows us to determine 
necessary properties for recovering the standard cosmological expansion in the late time limit. We have shown that 
the two-brane RS scenario cannot be made consistent with the desired expansion rate from standard 4D Einstein 
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gravity nor with observational constraints on Brans-Dicke gravity on the negative tension brane. This is true even if 
there is a stabiUzation mechanism, that is, constant or very slowly changing scalar field. A necessary condition for 
the recovery of a positive expansion rate is that the vacuum energy of the brane must be positive when the scalar 
field has been stabilized. 

The expansion rate of a cosmological solution can be changed in sign by introducing a conformal scaling of the 
metric. This is a highly contentious issue in the literature of scalar-tensor gravity theories. We have adopted the view 
that the induced metric must be that which governs the geometry of the spacetime (and hence the Hubble parameter 
of a cosmological solution). In section II D, we showed that, in any other frame, there would be second derivatives 
of the scalar field into the Einstein equation, resulting in an energy-momentum for the scalar field which violates 
positive definiteness. In addition, we see that, as with 4D Brans-Dicke gravity, the term in the Einstein equations 
proportional to the energy-momentum tensor still has a negative coefficient in the case of a negative tension brane, 
the possible change in sign of the Hubble parameter being due to the additional trace term. Dilatonic brane-world 
models are fundamentally different from 4D scalar-tensor theories in that the choice of frame is really an issue for 
the 5D theory, which then determines the frame of the 4D effective theory. In our calculation, we have started from 
a 5D Einstein frame and shown that the induced metric determines a frame which is preferable according to both of 
the criteria normally advanced in the study of 4D scalar-tensor theories, namely that matter follows geodesies of the 
metric and that the contribution to the Einstein equations from the scalar field is positively definite. 

Differences with results obtained by other authors who perform a dimensional reduction of the 5D spacetime by 
integrating over the bulk lead us to examine the difference between interpreting averaged or projected quantities as the 
relevant physical 4D entities. This is potentially an important philosophical issue in the study of higher dimensional 
spaces with warped geometries. We have argued that the averaging is difficult to define in a unique, coordinate-free 
way, whereas projection suffers from no such problem. In particular, it must be correct in the case that we have 
studied where the matter on the brane is localized by a (5-function, it may not necessarily be the case when the matter 
is defined in terms of excitations of higher dimensional gravity as in Kaluza-Klein models, for example. Clearly, this 
issue requires some understanding of how the matter is localized onto the brane from a fundamental point of view 
before it can be answered with any certainty. 

For a model based on a HW type dilaton, we have been able to obtain a next-to-leading order correction in (f) to 
the RS solution in the case of a universe dominated by radiation. The Friedmann equation is the same as for the 
RS model due to a cancellation of (j) dependent terms in the bulk contribution. Given that the HW type models 
have a rather different mechanism for achieving a cancellation of the effective cosmological constant, and that the 
asymptotic bulk geometry is likely to be somewhat different, it might seem surprising that the they reduce to the 
same effective Friedmann equation when the scalar field is assumed to be slowly varying. This suggests that form of 
the RS Friedmann equation is in fact universal in a brane-world which has any chance of being compatible with the 
standard Big-Bang model. 

It is natural to consider the case where the brane matter is non-minimally coupled to the scalar field. We discover 
that this does not change the cosmology when in a matter dominated era, but in the case of relativistic matter there 
is an additional term in the Friedmann equation. This gives rise to a cosmology where the universe expands from a 
finite value of the scale factor, with a period of accelerated expansion in the early universe. Although a toy model, it 
has the advantages of avoiding the initial singularity and flatness problem of the standard big bang model. Clearly, 
more should be directed in toward understanding this novel scenario. 
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